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Abstract 

We have implemented a systematic procedure for combining parton shower al- 
gorithms with next-to-leading order QCD calculations for the case of jet production 
in deep-inelastic electron-proton scattering. Using this method we have computed 
inclusive jet cross sections and jet shapes for the case of single-jet production and 
compared them to data from the ZEUS collaboration at HERA. We found good 
agreement between the data and our calculations, both for the jet shapes and the 
inclusive spectra. 
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1 Introduction 

X 

In the last two decades sophisticated event generators have been developed for all relevant 
physical scattering processes. For the case of eP scattering, generators like PYTHIA Q, 
HERWIG ||, LEPTO || or RAPGAP Q are available. In the simplest case the event generation 
procedure falls into three steps. First, a one-parton final state is generated from the 0(a Q s ) 
quark-parton-model (QPM) matrix element. Next, a parton shower (PS) algorithm §|| 
is attached to the one-parton final state. The PS takes into account arbitrarily high orders 
in a s , but only in the leading logQ 2 approximation as opposed to the exact treatment 
in fixed order matrix elements. This implies that the region of small k± around the 
original QPM quark is described quite well, which is important for the description of the 
internal properties of jets. Finally the PS is terminated at some lower cutoff Qq, where 
hadronization models such as the Lund String model || can be used to simulate the 
transition of partons to hadrons. 

The approximations involved in the construction of the PS allow to describe well only 
the region of small angle parton emission. The region of wide angle scattering should be 
more appropriately described by the hard QCD matrix elements of 0(a s ) which describe 
the QCD-Compton (QCDC) and boson-gluon fusion (BGF) processes with two final state 



partons. To include these processes into the event generators, the phase-space is separated 
into two complementary regions by means of some technical separation parameter, which 
we call Rtech in the following. In one ('soft') region, the QPM result plus PS is used 
to produce the cross section, whereas in the other ('hard') region the QCDC and BGF 
matrix elements are employed. The approximations involved in the PS approach suggest 
that the separation parameter R tec h should not be chosen too large in order to guarantee 
that wide angle scattering is appropriately described by the hard QCD matrix elements. 
One has to ensure that none of the partons generated by the PS populates the hard 
region above Rtech, since this would lead to double counting. Apart from separating 
the phase space into complementary regions, the parameter Rtech serves as a cutoff and 
ensures the fmiteness of the QCD matrix elements, which are divergent in the region of 
soft and collinear particle emissions. A problem occuring at this stage is that the overall 
normalization of the cross section is not yet fixed. For the special case we are discussing 
here, the total cross section can be calculated beforehand, which can be used to normalize 
the contributions coming from the QPM and QCD matrix elements. A lower boundary for 
Rtech is obtained by the requirement that the sum of the QCDC and BGF contributions 
does not exceed the total cross section. Otherwise Rtech can be chosen freely. 

The normalization of the individual contributions by calculating the total cross section 
beforehand is in general not feasible e.g. for the case of dijet production in eP scattering. 
It is therefore desirable to find a general, systematic approach to fix the cross section when 
fixed order matrix elements are combined with PS algorithms. To reliably calculate the 
normalization of the cross section, the complete next-to-leading order (NLO) corrections 
to the leading order (LO) processes have to be calculated. This includes the calculation of 
the real soft and collinear regions, as well as of the virtual corrections. The NLO correc- 
tions to the QPM graph in deep-inelastic scattering (DIS) have long been evaluated 0. 
Also for the DIS dijet case, several NLO calculations are avaible |8|,y,|ll],|ll|] , which employ 
either the subtraction method |8|, ^2|Jl3| or the phase-space slicing technique [p|,|Tl. |T5| . |T^] 



to treat the soft and collinear part of the real corrections. The problem of directly com- 
bining the NLO cross sections with the PS is the occurance of large positive and negative 
weights in the NLO calculation which makes it practically difficult to obtain numerically 
stable results. Furthermore, it is not straight-forward to use the hadronization models 
with negative weights. 

Recently, several proposals where made to combine NLO QCD calculations, including 



the virtual corrections, with PS's il7lJlq . |li4|20| . |2Tfl . In this paper, we will rely on the 



method proposed by one of us (B.P.) in ^D] and apply it to the case of single-jet inclusive 
cross sections in DIS eP scattering. The idea is to keep the separation of the phase-space 
into a soft and a hard region with help of the Rtech parameter, but to use the full NLO 
calculation in the soft region instead of the LO one. To ensure that the weights generated 
in the soft region are always positive, a method adapted from ]22] is employed where 
the s m j n parameter of the phase-space slicing method is chosen such that the sum of the 
Born, virtual and soft and collinear contributions vanish. In this way, the NLO corrections 
are calculated from the hard matrix elements, integrated within the soft region down to 
the s m in cutoff, which will always yield positive weights. While in the cutoff was 
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adjusted by hand, in [2(J it is calculated from the NLO corrections, thereby preserving 
the improved scale and scheme dependence of the NLO calculation. 

The outline of the paper is as follows. In section 2 we specify how to generate events 
with positive NLO weights and how we have implemented the method for the case of 
single-jet production in DIS eP scattering. In the following section we calculate inclusive 
jet cross sections and jet shapes as measured by the ZEUS collaboration at HERA for 
large Q 2 in the laboratory frame. We conclude with a short summary and an outlook. 



2 Event Generation with Positive NLO Weights 

In eP scattering 

e(k) + P(p) -> e(k') + X (1) 

the simplest hadronic final state consists of a single jet with a large transverse energy 
Et in the laboratory frame. The lowest order O{a s ) partonic contribution to this single- 
jet cross section arises from the QPM subprocess. At NLO, the single-jet cross section 
receives contributions from the real and the one-loop virtual corrections. The real cor- 
rections consist of the BGF and QCDC processes which are divergent for collinear and 
soft emissions. These divergencies are cancelled by the virtual corrections or are absorbed 
into the parton density functions (PDFs) of the proton. To enable a numerical treatment 
of the real corrections, one can employ the phase-space slicing technique and introduce 
a parameter s min which cuts out the singular regions. Within this method, the NLO 
cross section can be written as a sum of the one-parton final state up to 0(a s ) and the 
two-parton final state. The one-parton final state reads 

^haT^^mm) = <Tq G i J dx dPS {k ' +1) fi{x, fJi F ) (l + a a (fJL R ) JC q ^ q (s min , Q' 



=g,q 



M q -> q \ 2 . (2) 



The fi(x,fi F ) are the proton PDFs, the Lorentz-invariant phase space measure 
contains both the scattered electron and the partons from the photon-parton scattering 
process and the term o"o is defined as er = 4(na) 2 / \Q A xs). The are the charges of the 
quarks. The factor K q ^ q which depends both on the phase-space slicing parameter s min 
and on the hard scale of the process Q 2 contains the virtual and final state corrections 



and is specified in fl20fl . The function containing the initial state corrections is given 
by 

'Nr. 



C 4 MS (x,^,, mm ) = (^) 



Mx^ F ) \nl-^J+Br(x^ F ) 



(3) 



The functions A^x, fi F ) and B^^x, fi F ) are also specified in f20| . The final result, inde- 
pendent of s min , is obtained by adding to the one-parton contribution the contribution 
containing the two-parton final state, integrated over those phase-space regions where any 
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pair of partons i,j has Sy > s m i n with sy = (pi + pj) 2 : 



a 



2parton / \ 
had l s minj 



E e2 i J dx ^PS (fc ' +2) 4vra s (/i i? ) 



fi(x,fi F )] \M ( 



'JK°min 



+ ifg(,X,fi F ) \ M g^qq\ 2 



The condition 



j lparton 

had f^?J = 



dx dQ 
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can be used to determine the function 20 



s mia(^J X i Q A ) = ex P 
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and 



(9) 
(10) 

(11) 



i=q,q 



Inserting the s^° n function (||) into Eqn. (f|) as a lower integration boundary for each 
phase-space point (x, Q 2 ) will give the complete answer for the total cross section at 
NLO. It is important to note that the s£j° n function depends on the factorization and 
renormalization scales, so that the improved scale dependence of the NLO cross section 
as opposed to the LO one is preserved. As was studied in detail in [PD| , the s^° n function 
Eqn. (^j is small enough for the soft and collinear approximations to remain valid which 
are made to evaluate the 0(a s ) one-parton final states. The NLO single-jet cross section 
calculated within the standard approach could be reproduced within 1-2 % with this 



method [20 



For the calculation of the cross sections in the next section, we have implemented 
the function s^° n into the DISENT program |J and combined the NLO cross section thus 
obtained with the PS algorithm from PYHTIA [l]. The reason we chose DISENT is to be 
able to easily extend the results obtained here to the dijet case later on. In the dijet 
case also the subtraction terms from the subtraction method have to be used (see f2(| for 
details). The steps we have performed to generate events are the following: 
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• Define the soft (PS) region by < R tec hQ 2 for at least one pair of partons 

• The hard region given by Sij > RtechQ 2 for all partons i,j is described by the hard 
two-parton matrix elements (BGF and QCDC). We do not attach the PS to the 
partons of the hard region, although this could in principle be done. 

• For each (x, Q 2 ) calculate s^? n according to Eqn. @. When < s£]? n for any one 
pair of partons reject the event. When > s£j° n for all start the PS 
algorithm with the four- vector of the QPM quark as input. The weight associated 
with this event is the one obtained from the two-parton final state. These weights 
are positive by construction. 

• Reject any partons from the PS that lie outside the PS region. 

The parameter Rtech should be larger than s^f n /Q 2 but not too large to ensure that 
the hard region is appropriately described by the fixed order matrix elements. Rtech is 
typically of the order of 1. We have checked that our cross sections are not sensitive to 
the exact value of Rtech by varying it by a factor of 2. The starting scale for the PS is set 
by the matrix element cutoff ||. 

We have not yet implemented the initital state PS, but only the final state PS. The 
parton emission due to the initial state PS will populate mainly the forward region, which 
is excluded in the data with which we are comparing in the next section through kine- 
matical cuts. Therefore we believe that the effects of the initial state PS are small. The 
evolution of the incoming parton is taken into account through the evolution equations for 
the partons in the proton. Furthermore, hadronization effects have not been taken into ac- 
count in our calculations. It is however feasible in our model to implement hadronization, 
since it can be naturally attached to the PS. 

For all calculations in the next section we have set the renormalization and factoriza- 
tion scales equal to Q 2 and employed the CTEQ4M PDFs for the proton. 



3 Results 

The ZEUS collaboration has measured inclusive jet cross sections in neutral current (NC) 
DIS eP scattering at y/s ~ 300 GeV for Q 2 > 125 GeV 2 f23p|. The phase-space of 



the electron has been further reduced by restricting its energy E e > and the inelasticity y e 
to E e i > 10 GeV and y e < 0.95, respectively. The jets were reconstructed using a kx 



cluster algorithm [p5j in the laboratory frame. The reconstructed jets were required to 
have a minimum transverse energy of > 14 GeV and a pseudorapidity in the range 
— 1 < rj < 2. The data are corrected for detector effects. 

In Fig. [I] we show the comparison of our calculation, given by the full line labeled 
'DISSET', to the ZEUS data in four regions of Q 2 , namely Q 2 > 125,500,1000 and 
2000 GeV 2 . We see an overall good agreement. Since our predictions are on parton 
level, we have checked the influence of the hadronization using the LEPT0 event generator 
@. We found that the overall effect of the hadronization was small, below 5%, with a 
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tendency to lower the parton level predictions. We have furthermore checked that the PS 



changes the distribution only marginally, also below 5%. In [24 the data were com- 
pared to a standard NLO calculation with DISENT 0. We have redone these calculations 
for comparison and also plotted the results in Fig. [I], shown as the dashed line labeled 



'standard DISENT'. Our results are in agreement with those from |24| . The agreement 
between data and NLO theory is similarly good as for our calculation. From this compar- 
ison we deduce that our calculation reproduces correctly the standard NLO result, which 



confirms the findings in ||20|| . We emphasize that the correct normalization of the cross 
sections is a non-trivial result of our procedure. We obtained the normalization without 
using information from the total eP scattering cross section. The normalization comes 
directly from the NLO matrix elements, modified by the sj^° n function. 

A similarly good description of the data is seen in Fig. ^| for the rj distribution, which 
is integrated over E? > 14 GeV and Q 2 > 125 GeV 2 . The data are described rather well, 
both with the standard NLO calculation as well as with our new calculation. 

We proceed to a comparison of our calculation with jet shape measurements. The 
ZEUS collaboration has measured the differential and integrated shape of jets in neutral 
current DIS events with Q 2 > 100 GeV 2 for jets with E? above 14 GeV and — 1 < rj < 2 
The jets are reconstructed using an iterative cone algorithm in the (rj, <p) plane 



27| , |28| . For the scattered electron further restrictions are E e * > 10 GeV and y e < 0.95, 
as for the previously discussed data. The differential jet shape is defined as the average 
fraction of the jet's transverse energy that lies inside an annulus in the (rj, 0) plane of 
inner (outer) radius r — Ar/2 (r + Ar/2) concentric with the jet defining cone p9| |: 



^^^ ^It)^' (12) 

l\ je t s Ar jets E T [[) } H) 

Here, E?{r — Ar/2, r + Ar/2) is the transverse energy within the given annulus and Nj ets 
is the total number of jets in the sample. The differential jet shape has been measured 
for r values varying from 0.05 to 0.95 in Ar = 0.1 increments. The integrated jet shape 
is defined by 

Wjets jets E T {U,K) 

By definition, ^(R) = 1. It has been measured for r values varying from 0.1 to 1.0 also 
in increments of Ar = 0.1 

In Fig. |]we compare our calculation (full line) and a standard NLO calculation (dashed 
line) to the ZEUS data which are corrected to hadron level. The description of the jet 
shape within our approach is rather good. The jets produced are slightly too narrow, our 
calulation being about 8% above the data in the lowest r bin. We have checked with LEPT0 
that the hadronization effects can account for the difference between data and prediction, 
i.e. including hadronization in our calculation would bring our results to good agreement 
with the data. We furthermore see that the standard NLO prediction is much too large 
for all r. This implies that the jets produced by the standard calculation are not broad 
enough. Reducing the cross sections by taking into account hadronization is not sufficient 
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for the standard NLO calculation. This result is not surprising since the standard NLO 
calculation provides only a LO prediction for the jet shape in which at maximum two 
partons are combined to give a jet. 

In Fig. |] we finally compare our calculations to the ZEUS differential jet shapes for 
four different Et regions. We find similar results as for the integrated jet shapes. The 
standard NLO calculation (dashed line) is clearly too narrow in all four Et regions, 
whereas our calculation (full line) gives a much better description. In the lowest Et bin 
our calculation produces slightly too narrow jets, whereas in the largest Et bin the jets 
tend to be too broad. Using LEPTO, we have checked the hadronization corrections also 
for the differential jet shapes and found the discrepancies between our calculation on the 
parton level and the data can be accounted for by hadronization effects. 



4 Summary and Outlook 



We have investigated the method described in |20| to combine a fixed NLO QCD cal- 
culation, including the real soft and collinear as well as the virtual corrections, with a 
leading log PS. We have implemented the method in DISENT and attached the PS from 
PYTHIA. We have compared our calculation with data from ZEUS, both for inclusive jet 
cross sections in the laboratory frame, as well as for jet shape measurements. For both 
measurements, good agreement with our model was found. In addition, good agreement 
was found between the standard NLO calculation and our results for the case of the inclu- 
sive spectra in transverse energy and rapidity, where the influence of the PS is marginal. 
This shows that we have obtained the correct NLO normalization of the cross sections 
together with a good description of the internal structure of the jets. 

To obtain a full event generator, the initial state PS as well as hadronization need to 
be implemented in our program package. This will be done in the near future. The next, 
more complicated step is to proceed to the two-jet final states and combine the NLO 
calculations for this case with the PS. Although some details have to be clarified for this 
case, we could show the principle feasibility of our method. 
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Figure 1: Differential cross sections dcr/dE? for inclusive jet production in NC DIS events 
integrated over -1 < r] < 2 for Q 2 > 125, 500, 1000 and 2000 GeV 2 . The ZEUS data are 
compared to our calculation labeled 'DISSET' (full line) and to a NLO calculation, labeled 
'standard DISENT' (dashed line). 
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Figure 2: Differential cross sections da/drj for inclusive jet production in NC DIS events 
integrated over E T > 14 GeV and Q 2 > 125 GeV 2 . The ZEUS data are compared to our 
calculation (full line) and to a NLO calculation (dashed line). 
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Figure 3: Integrated jet shape \I/(r) in NC DIS events integrated over Q 2 > 100 GeV 2 for 
jets with E T > 14 GeV and — 1 < rj < 2. The ZEUS data are compared to our calculation 
(full line) and to a NLO calculation (dashed line). 
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Figure 4: Integrated jet shape p(r) in NC DIS events integrated over Q 2 > 100 GeV 2 for 
jets with — 1 < i] < 2 in four different ranges of Et- The ZEUS data are compared to our 
calculation (full line) and to a NLO calculation (dashed line). 
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